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Introduction 

Tine purpose of this paper is to develop and to illustrate by a number 
of exaiiples some interesting and useful extensions of a procedure, gener- 
ally known as Dunkerley's method, for estimating the gravest frequency of 
an undamped linear mechanical vibrating system. The theory is developed 
in the next section hereof. The equivalent of equation 15 was given in 
the year 189^ by S. Dunkerley (3) who obtained the result enpirically 
based on calculations performed while investigating the vibrations of 

shaft and disk systems; he certainly did not recognize the generality or 
the theoretical basis of his formula. Tenple and Bickley (12) discussed 

the procedure in 1933 indicating its applicability to both limped and 
continuous mass distributions. However Bickley and Talbot (1) In a 
later (1961) textbook on vibrations do not mention the method. South- 
well (11) treats the method in his well known treatise (1936, 19^1) without 
citing Dunkerely. This may be the reason that the method is sometimes 
called the Dunkerley-Southwell method. However, although Southwell may 
have arrived at the result independently (he did indeed discover another 
method of determining lower bounds; cf. Lamb and Southwell (6)), .his book 
does cite Temple and Bickley (12) among the general references. 



- 2 - 



A formulation in terms of integral equations and their eigenvalues, 
iterated kernels, etc., which is appropriate for continuous mass distrib- 
utions, is a logical extrapolation, by analogy, from the formulation in 
terms of matrices, their eigenvalues, powers of matrices, etc., which is 
appropriate for systems having a finite number of degrees of freedom, so 
that, in a sense, all the work reported herein stems from Eunkerley (3) 
and Temple and Bickley (12). Indeed it was via this path that the writer 
was led to the results reported here. However, the relationship with the 
theory of integral equations, a thoroughly developed and explored discipline, 
is so very close that it seemed likely that the developments at which the 
writer has arrived had been anticipated, in purely mathematical context, by 
an earlier writer. Indeed this is the case, Mikhlin (8) definitely attri- 
butes the ideas to Mikhlin (7). In neither of these references, however, 
does Mikhlin refer to EUnkerley (3), Temple and Bickley (12), or South- 
well (11). It seems clear that Mikhlin arrived at his results without 
being aware of Dunkerley's formula or of its relation to his own work. 

Accordingly, the writer believes that it is appropriate to call the ex- 
tended procedure by the name Dunkerley-Mikhlin. The present paper may 

be considered as a brief exposition of this method, and some extensions, 
and of its application to a variety of problems of engineering vibration. 

Derivations 

We consider an undamped linear vibrating system characterized by 
N*N symmetrical matrices K (stiffness) and M (mass). Suppose that there 
are p rigid body modes which are known (by inspection). (Thus, if p > 0, 

K has no inverse.) There exists a modal matrix U, not necessarily unique,' 
and a (diagonal) spectral matrix ft 2 in which frequency-squares are arranged 
in non-descending order, such that 
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(la,b) 



KU » KUfi 2 ; U T MU = I = N*N unit matrix 
Now consider any other system which is similar to the original system 
except that p additional constraints have been incorporated so as to eliminate 
the rigid body modes. Let the flexibility of this system be C. 

' Consider the kth mode (k>p) of the original system. There are constants 
a. , such that 

jk P 

\ ■ \ + 2>jk u j (2) 

does not involve applying loads to the added constraints. Thus, for k>p, 



v k - CKv k - 



P 

:£ 



au k + c jti a jk Ku j " = “w 



(3) 



Define the filtering matrix 



I - Su.uTm 

0-1 J J 



W 



which is such that 



% 



0 If kip 
u^ if k>p 



Thus, for k>p, 



■ ““v = “k \ 3 “k \ + “k £ a jk u j 

j=i 

while for k^p, 

cmf^ = 0 

Equations 6 may be combined to give 
UA + B ■ CMFU 

where 



A = diag[0 0 ... 0 wi 



p+1 w p+2 *•* 






(5) 



(6a) 



(6b) 

(7) 

( 8 ) 
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and B is a matrix the first p columns of which are zero and the remaining 



colunns of which are linear combinations of Uj, i^, 
the matrix iT^B = U^MB, the kth column of which is 





(9) 



and, because of equation lb, this can have nonzero elements only in the 
first p positions. Thus A + U^B and (A + U -1 B) n are of the forms 




respectively, where 0, is a p*p matrix of zeros, 0 2 is a (N-p)xp' matrix 
of zeros, E is the nonzero part of U _1 B, and D is the nonzero part of A. 
Finally, from the equality 



by taking traces, we get 
N 

y> 2n = tr[(A+ ir'B) n ] = tr[(U-'CMFU) n ] = tr[(CMF) n ] = 
k=p+l “ 

where Q is defined as the matrix triple product CMF. Note, however, for 
computational convenience, that 



We have used the fact that if A and B are matrices conformable in either order 
then tr(AB) = tr(BA). 

V/e will refer to equation 12 as the nth DM (Dunkerley-Mikhlin) 
evaluation and to the corresponding relation 



(A + lT l B) n = (U^CMFU) 11 



( 11 ) 



tricf 1 ) = tr[(CMF) n ] = tr[(MFC) n ] = tr[(FCM) n ] 



(13) 



Vi^ Ctr(Qn ) r 1/2n d^) 

as the nth DM approximation. We will repeatedly ustthe symbol § to mean 
"greater than but approximately equal to," and a symbol indicating the 
opposite order. 

If M is diagonal, if there are no rigid body modes (i.e,, if p=0) , 
and if we take n=*l, equation 14 becomes ^ 

«i* » tr(Q) = tr(CM) = (15) 

This is what is generally referred to as Dunkerley's formula. It was 
(obliquely) stated in 189^ by S. Dunkerley who regarded it as an empirical 
representation of calculations he had made of shaft and rotor frequencies. 

Simple extensions to continuous mass distributions are obvidus and 
have been discussed by many writers, among the first of whom were Temple 
and Bickley, Reference 6. However, the theoretical basis for the contin- 
uous case seems to have been first established by MLkhlin, who, however, 
did not consider problems of mechanical vibration and who seems to have been 
unaware of IXmkerley's work. In the continuous case, matrix multiplication 
is replaced by integration and the compliance matrix C is replaced by a 
symmetrical function of two variables, z(x,y) which gives deflection at x 

due to unit loading at y. This notation is appropriate to a one-dimensional 

herein 

field; however, examples given irv Reference -1- illustrate cases involving 

two dimensional fields. The mass matrix M which heretofore need not have 

been a diagonal matrix, now becomes the equivalent of an infinite diagonal 

matrix, namely, a function m(x) specifying mass per unit length. The matrix 

Q is replaced by the function p L 

q(x,y) ** m(y)[z(x,y) - 2»u. (y) f u.(0 m(£) z(x,0 d£] (16) 

k=l K *0 K 
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where the functions u^x), k=l,2,. . . ,p, describe the known orthononnal 

rigid body modes, satisfying 

L 

| m(0 Uj(0 dZ = 6^ (17) 

0 

(Kronecker delta). In equation 1 6, z(x,y) is the conpliance of any 
system which is like the given system but with additional constraints 
so as to eliminate the rigid body modes. The function 

z(x,y) is sometimes called the Green's function. 



The continuous analog of raising the matrix Q to the nth power 
is the formation of the nth iterated function 



q^y) = 



q n _ 1 (x,0 q(£,y) d£; q 1 (x,y) = q(x,y) 



and the operation corresponding to taking the trace of Q n is 

co 

t L 



k=p+l • o 



k=p+l 

If there is both lumped and distributed mass, we may write 

r 



m(x) + 



m(x) = m(x) + L-j m 5(x-x, ) 
k-1 K K 



( 18 ) 



(19) 



( 20 ) 



where 6(x-x^) denotes the Dirac "function," m(x) denotes a continuous mass 
distribution, and the sum represents r distinct point masses m^ located 
at x = x^, k=l,2,.,.,r. 



Illustration 

We consider a finite element model of a segment of an Euler-Bemoulli 

T t 

beam. The displacement vector is u * [x, x 2 x 3 x 4 ] where x t is the (up- 
ward lateral) displacement at the left end, x 3 is the similar displacement 
at the right end, x 2 /L is the slope at the left, and x 4 /L is the slope at 
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the right. The length of the segment is L. The FEM consistent stiffness 
and mass matrices are 



K = (2EI/L 3 ) 


'6 3-6 3 


; M » (m/420) 


156 22 54 -13 




3 2-31 




22 4 13 -3 




-6 -3 6 -3 




54 13 156 -22 




.3 1-3 2 




-13 -3 -22 4. 



First we exhibit the classical solution. The modal matrix (the first two 



colums of which are not uniquely determined) is 



U 



lab c 

0 «*2a -6b —12c 

1 -a b -c 
0 -2a 6b -12c 



/\/m t where a 2 = 3, b 2 = 5, and c 2 = 7. 



and the spectral matrix is 

n 2 - diagCw 2 w 2 G) 2 w 2 ] - (840EI/mL 3 ) diag[0 0 6/7 10] 

2 2 3 4 

There is no difficulty in verifying equations la and lb. 

Next, we enploy the EM method and verify the correctness of 
equation 12. Vfe constrain the system by rigidly fixing the left end 
so as to form a cantilever. We easily determine 



C - (L 3 /6EI) 


'0 


0 


0 


o' 


; CM » (mL 3 /2520EI) 


' 0 


0 


0 


o' 




0 


0 


0 


0 




0 


0 


0 


0 




0 


0 


2 


3 




69 


17 


246 


-32 




0 


0 


3 


6 




.84 


21 


336 


-42. 



and we also calculate 



F « 



-3 -4 



3 



1 



/30 



36 33 36 3 
3-1 -3 4 
.36 3 -36 33. 



Next, we obtain 



Q * CMF 



(mL 3 /25200EI) 



' 0 0 
0 0 



-3 -19 



0 -35 



0 O' 
0 0 
3 16 , 
0 35 
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Finally we calculate 



tr(Q) = 38mL 3 /25200EI = ui" 2 + w 4 ' 2 
tr(Q 2 ) = 12 34 ( rnL 3 /25200EI ) 2 = to" 4 + u>‘ 4 
tr(Q 3 ) 3 42902(mL 3 /25200EI) 3 3 ui' 6 + o^ 6 
and so on. All of these calculations check out correctly. 

One may observe that the two diagonal elements of Q give directly 
w‘ 2 = 35mL 3 /25200EI , oj~ 2= 3mL 3 /25200EI 

but this Is only fortuitous In the present illustration as may be seen In 
examples lb and lc which appear later. 

Now we illustrate the equations employed in the derivation of 
equation 12, First, it is easy to verify that 

FU 3 [0 0 u 3 u^ ] 

Next we form 



0 0 

-1 -1 

0 -1 



0 

0 

1 

0 



0 “ 

0 

0 

1 



and obtain 

[v 3 V 4 ] ■ CK[u 3 U 4 ] 



0 0 
0 0 
6b 10c 
12b 0. 



/Jm. 



which obviously satisfy the constraints which have been introduced at the 
left end. We can also calculate the coefficients a 13 = 2b, a 23 * -3b/a, 



l l4 



5c, and a 3 -6c/a, so that, indeed 



2 4 

v 3 - u 3 + a 13 u 3 + a 2j u 2 ; v 4 = u + a. u. + a. . u. 



4 14 1 24 2 

We can also calculate and verify that 

WjCMu, = ( 6/7 ) ( 84 OEI/mL 3 ) ( bmL 3 /2520EI ) [ 0 0 21 42] T /^m = v 3 



with a similar calculation relating u 4 and v 4 . 
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Next we calculate 



CMFU » QU * (v / rnL 3 /840EI) |"0 0 0 0 

0 0 0 0 
0 0 7b c 
.0 0 14b 0 

A - (mL 3 /25200EI) diag[0 0 35 3] 

UA « (VmL 3 /25200EI) fo 0 35b 3c 

0 0 -210b -36c 
0 0 35b -3c 

_0 0 210b -36c 

B - CMFU - 



UA - (»^nL 3 /25200EI). 


0 


0 


-35b 


-3c 




0 


0 


210b 


36c 




0 


0 


175b 


33c 




0 


0 


210b 


36c 


- (mL V25200EI) [0 


0 


35(a u 


, U l +a 2 



Noting that U" 1 
U* l QU = 



T 

(MU) , we finally calculate 

(mL 3 /25200EI)f0 0 70b 15c 

0 0 -3.5ab -6ac 

0 0 35 0 

0 0 0 3 . 



A + U" l B 



so that 



D = (mL 3 /2520QEI) 


‘35 o' 


; E = (mL 3 /25200EI) 


70b 15c ' 




: 0 3 




_-3.5ab -6ac. 



We also verify that 

tr(U _1 QU) - 38mL 3 /25200EI, tr[(U- l QU) 2 ] » 1234(mL 3 /2520QEI) 2 , etc. 
Although It is worth remarking that neither K nor C possesses an 
Inverse, this is to be expected In cases for which there are rigid body 
modes | and the modified system involves perfectly rigid constraints. It 
is of some interest to note the form of the product CK. 
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Applications 



Example l t part a. A uniform, massless cantilever beam has concentrated 
masses m, 9m, and 4m at distances 9L, 21 L, and 27L, respectively, from 



the fixed end. It is desired to estimate the lowest frequency of harmonic 
oscillation. Taking the lateral deflections of the point masses, namely 



x^, x 2 » and x^, as the elements of a vector u, we easily determine 
M = m dlagCl 94]; 



C 



(9L 3 /EI) 



27 81 108' 
8l 34 3 490 
108 490 729 



Q » CM = (9mL 3 /EI) 



■ 27 729 432 ' 
81 3087 I960 
108 4410 2916 



Prom equation 12, taking n=l, we get 

V V 

(EI/54270mL 3 ) 2 = 4.2926*10" 3 (EI/mL 3 ) 2 



Taking n=2, v/e get 

(E 2 I 2 /2878o89084m 2 L 6 )^» 4.3174l6'10" 3 (EI/mL 3 ) i ' 2 
A Rayleigh approximation, details of which are not given here, gives 
the numerical coefficient for an upper bound as 4. 317542 *10“ 3 , so that 
the first approximation is in error by less than 0.6# and the second 



is in error by less than 0.003#. 



Example 1, part b. Now suppose that the fixed support at the left of 



this beam is replaced by a frictionless pivot. The frequency co l = 0 
corresponds to the rigid body mode u, s [3 7 9]^/(774) /2 . We obtain 



F - 



765 -189 -108' 
-21 333 -252 
-27 -567 450 



/774; 



Q 



= CMF = (9mL 3 /86EI) 



’ -702 -810 864* 
-6198 -10962 10952 
-9858 -22806 21024 



and we have no difficulty in obtaining 

w a £ 0.03195(EI/mL 3 )^ 2 , 0.033494(EI/mL 3 )^ 

for the first and second DM approximations, respectively. A Rayleigh 
approximation gives the numerical coefficient as 0.033586 so that the 
errors are less than 5# and 0,3# respectively. 



Example l t part c. Now suppose that the beam is completely free at the left 
so that there are two rigid body modes, involving translation and rotation. 

We use the same matrix C as above although a simpler one could be constructed 
since the leftmost segment of length 9L obviously does not enter the present 

ip y m 

problem. Vfe take u t - [1 1 1] /(l^) and as a combination of [1 1 1] 

T 

and [3 7 9] proportioned for orthonormality. Thus we find Uj = . 

[15 1 -63 T /(378)^. Vfe calculate 



F - 


‘ 42 -126 84' 


/126; Q - (96mL 3 /EI) 


'0 0 O' 




-14 42 -28 




1-3 2 




21 -63 42 




.3 -9 6 



and we get tr(Q) » 288mL 3 /EI, tr(Q 2 ) ® 82944(mL 3 /EI ) 2 both of which give the 
exactly correct result 

» (EL/288mL 3 ) ^ 

since there are only three frequencies the smallest two of which are known 
to be zero. 

Example 2. Now we consider a cantilever beam in the form of a truncated 
right circular cone having base radius b, tip radius a, and length L, and 
having uniform physical properties. Assuming an Euler-Bemoulli model, 
we obtain the bending compliance 

z(x,y) 3 US-n)[(l-an) -2 - (l+2an)] + 2a 2 n 3 /(l-an))/6Ba 2 
where 

B = TrEbV4 = ^ EI ^root* a “ ( b “ a ^/ bL J £ 3 Max(x,y}; n = Min(x,y} 

The mass per unit length is 

m(x) » mo(l-ax) 2 ; mo = mb 2 y 

mo being the value pertaining to the root end and y denoting mass density. 

Performing the indicated evaluations gives 

; L 

w* 2 » j m(x) z(x,x) dx = moL 1 * (l+4q)/60B 

Jo 

for the first estimate and, noting the symmetry of z(x,y), 
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- rno 2 L 8 [3(q 4 ~l)/8 - q 4 log e (q)/4 + 3p 2 - 5p 3 + llp 4 /^ 



- 3p s /5 - p 6 /20 - p 7 /70 + 9pV80 - p 9 /5 + 2p lo /25]/l8B 2 p 0 



where q = a/b, p = 1-q. Although both of these evaluations give results for m, 
which are smaller than the correct value, the second gives a result which 
has a maximum error of 0.8# (for q = 0) for 0 = q = 1. The first, and 
much simpler, evaluation is much less accurate, but it serves as a guide to 
the following approximation 



which has a maximum error of for q = 0.1. The "exact" results, for 
comparison, were calculated from formulas given by Conway, Becker, and 
Dubil, Reference 2, 

Next, we consider some cases for which the mass distribution is 
of the form described by equation 20. The first and second DM evaluations, 
respectively, take the forms 



The first term in these evaluations is what we had previously for continuous 
distributions and the last term may be identified with the trace of the 
appropriate matrix. 

Example 3. Suppose that the cantilever beam of example la itself has a 



Wj 8.72[(l~0.0l6q)/(l+5.053q)] J/2 (B/m o L 4 ) ,/2 




r 
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mass 4m distributed uniformly along Its length. (Note that the upper 
limit of integration is 27L.) We havp m(x) = 4 jtv' 27L and z(x,y) = 
(3xyw-w 3 )/6EI, where w = Min{x,yj. We calculate 

w. £ [(W27L)(27L)Vl2EI + 54270mL 3 /EI] = 0.004054(EI/mL 3 ) ^ 

* 

for the first DM approximation. Note that we have seen the number 
54270 in example la. The second EM approximation is more labor; it 
gives 

w, £ [ll«3 ,7 /35 + 4759166988/7 + 2 878089084 ]"^(EI/mL 3 ) 1/$ 

= (81*1554923033/35)' ^(EI/mL 3 )* 2 = 0.00408289(EI/mL 3 )^ 

An "exact” evaluation by transfer matrix procedure gives the numerical 
coefficient 0.00408305, so that th6 errors for the first and second 
approximations are 0.7$ and 0.004 % respectively. 

Exarple 4, Consider the case of a uniform cantilever beam of length L 
and mass pm having a concentrated mass (l-p)m at its tip. The compliance 
function z(x,y) is given in sample 3 and the mass function is m(x) = 
pm/L + (1-p) m <$(x-L), The first and second DM approximations give 
[(3EI/mL 3 )/(l-3p/4)]^; {(a):-1.5$ at p - 1} 
o)j> [5040(EI/mL 3 ) 2 /(560-856p+329p 2 )] 1/4 ; {(b):-0.017* at p - 1) 

( In these and other estimates given in this example 4, the estimate is 
distinguished by a lower case letter identification, followed by the maximum 
percent error and the value of p for which it obtains.) 

This example affords an opportunity to discuss other types of simple 
but accurate approximations. A Rayleigh approximation based upon the 
deflection function y » 3Lx 2 -x 3 gives 

Wj £ [(3EI/mL 3 )/(l-107p/l40)]^; {(c):+1.5$ at p - 1) 

This is an upper estimate while (a) is a lower estimate. An obvious average 
is 

«*>i ~ [( 3EI/mL 3 )/(l-53p/70) ] ; {(d):-0.39$ at p - 0 . 83 } 
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Approximation (d) Rives equal weighting to (a) and (c). This happens 
to be approximately optimal in this case but this is only fortuitous. The 
coefficient of p should lie somewhere between 105/1*10 and 107/1*10. Both are 
exact for p - 0. We can choose the coefficient so as to make the result 
exact for another value of p t say p = 1. Ibis is a uniform cantilever for 
which 

to, = (1. 87510*10687... ) 2 (EI/mL 3 )^ 

Thus we obtain 

Uj » [(3EI/mL 3 )/(l-0.757328p)] K2 ; {(e):-. 366% at p = 0.81} 

Exactly the same result is obtained by replacing the first term in the 
first DM evaluation by the known exact value appropriate for a uniform 
cantilever with no added mass at the tip. This suggests replacing the 
coefficient 11/1680 in the first term of the second DM evaluation by 
(1.8751.. .) -8 and we get 

a, « [630(EI/mL 3 ) 2 /(70-107p+*11.1222805p 2 )] 1/< ‘; { (f) :-0.006$$ at p = 0.86} 

This accuracy is almost incredible. It goes without saying that the 
accuracy of these approximations is far greater than that of the physical 
theory to which they pertain. 

Example 5. \js next take up an example which combines the features of 
mixed (i.e., lumped and distributed) mass distribution and rigid body 
motion. We consider the axial motion of a uniform elastic bar of mass m, 
length L, and axial stiffhess EA, which has a concentrated mass 2m at the 
left end and a concentrated mass 3m at the right end. This dumb-bell 
shaped object is not constrained or tied down at any point. The mass 
function is 

m(x) = m/L + 2m6(x-0) + 3m6(x-L) 

We consider the compliance of a system fixed at the left end. 

z(x,y) = Min(x,y)/AE 

-15- 



and we obtain 

q(x,y) = [m(y)/AE][Min(x,y) - ( 8 Lx-x 2 )/12L] 

Thus the first DM evaluation gives 

w, > (9AE/13mL)^ = 0.832(AE/mL)^ 
and the second gives 

w >(l44A 2 E 2 /242m 2 L 2 )^* = 0.87829(AE/niL) /,/2 
The exactly correct coefficient is the smallest positive root of the 
equation tan x = 5x/(6x 2 -l) t which is approximately 0.87935. Thus the 
errors are 5 .^% and 0 . 12 $ respectively. 

Exanple 6 We determine the fundamental 
harmonic frequency of a uniform hinged- 
guided beam of length L. The mass func- 
tion is m(x) = mo » const. The conpliance 
function is 

z(x,y) = z(y,x) = [3xy(2L-y) - x 3 + <x-y> J ]/ 6 EI 



r 



Sketch of uniform 
hlngea-guided beam 



where 



f (x-y ) 3 if x=y 
<x-y > 3 » ■{ 

[0 if x*y 



The DME1, given by equation 19 with n 3 1, is 

r L 



to 



f £ SjW* 2 = | moZ(x,x)dx = moLV6EI; 2.^95(EI/m 0 L 4 /* 



The DME2 is obtained from equation 19 with n 3 2, viz. 

wT^ « = [ f q(x,C)q(S,x)d£dx = 

i=l 1 > 0*0 

» 2 [ [ m(x)m( 5 )z 2 (x, 5 )dCdx 

> 0*0 

because of the symmetry of z(x,y). The important consequence of this 

is that we can take 

z 2 (x,0 - [3xC(2L-x)-^ 3 ] 2 /36E 2 I 2 

and not have to concern ourselves with the alternate form applicable 
for £ > x. Continuing the evaluation we arrive at ( 17 / 630 ) (moL 1 * /El) 2 , 
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from which we obtain 

w x £ 2. i J673(EI/m 0 L") V51 
Ihe exact result Is 

u x = (Tr 2 /4)(EI/m 0 L‘*) Vi = 2.^674(EI/m 0 L'*) ,/X 
so that the excellence of the approximations is evident. 

Having the excellent approximation to w x given by the DMA2, we 
can get a reasonable approximation to the second frequency, viz. 

wl 2 c: - C(“T [, ) EMA 2 ^ !B [1/6 “ (17/630)' /2 ]m o L ,, /EI 

giving 

o) 2 ss 20 . 4 (EI/moL 1 * /* 

Ihe correct coefficient is (3V2) 2 = 22.2. 

This example permits us also to exhibit another possibly useful 

feature of the method. In this case it is known that 

= [(k-£)*] 2 (EI/moL 1 * )*\ k=l,2,..., 

-2 -4 

so that our precise evaluations of Eto^ and Ew^ lead to 

oo co 

^~ l . n"^ 1 = irV9 6; ! n”^ » 17fr 8 /l6l280 

n=l,3,5,.. . n=l,3,5>... 

These are well known results; cf. Jolley (5). However, for other cases 
the corresponding sums may be unknown. For this reason we have looked 
into all possible cases of a uniform beam having end conditions of the 
following types: fixed, hinged, guided (see right end of the figure), and 
free. Letting 

n 1 * - w 2 m 0 L , ’/EI 

the results may be exhibited in the form shown in Table 1. 

In Table 1 the physical case described by the specification of end 
conditions corresponds to the frequency equation indicated. Ihe sum of 
reciprocal second and fourth powers of the frequencies (except for factors 
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Table 1 Some results for uniform beams 



CASE 


EQUATION 


END CONDITIONS 




Z u” 1 * 




1 


cosn 3 sechn 


clanped-clanped ( * ) 
free- free 


1/420 


71/17463600 


4.7300407449 


2 


tann 3 tanhn 


clanped-hingsd( *) 
hinged-free 


1/210 


13/727650 


3.9266023120 


3 


tanq = -tanhn 


clanped-guided ( * ) 
guided- free 


1/30 


29/2835 


2.3650203724 


4 


cosn 3 -sechn 


clamped- free 


1/12 


11/1680 


1.8751040687 


5 


sinn 3 0 


hinged-hinged ( * ) 
guided-guided ( * ) 


1/90 


1/9450 


3.1415926536 


6 


cosn 3 0 


hinged-guided 


1/6 


17/630 


1.5707963268 



* The root, n 3 0, is to be excluded for cases marked with an asterisk 



of EI/n^L* and E 2 I 2 /m 0 2 L 8 , respecively) are shown. The smallest root 
of the frequency equation is also given. This information may be useful 
in augnenting infomation contained in work by Young et al. (16) and by 
Gorman (4), It may be remarked that neither of these references treats 
the end condition ^ = 0 which we call "guided;" this condi- 

tion is certainly not technically important but its inclusion seems to 
be indicated on the grounds of completeness. 

In the table, an asterisk indicates that the obvious root, n = 0, 

Of the frequency equation is to be excluded. The sums include only non- 
zero frequencies, and indicates the smallest positive root. We 
evidently have such results as 

71/17463600 

(i) imi 

where denotes the ith positive root of equation j. From this 
evaluation we could obtain 

n ^4. 7191 
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Exanple 7 In order to illustrate a case for which the field of integration 
is two dimensional, we consider a uniform circular elastic plate which is 
clanped at its outer radius r = a. The influence relationship, z(P,Q), 
is given by Timoshenko (1*1) in the form of a series, but for our purposes 
all that is needeiis the "self" influence, i.e., the deflection under the 
unit lateral load 

z(P,P) = (a 2 -r 2 ) 2 /l67rDa 2 

which is also given in (1*0. Thus, with mo denoting mass per unit area of 
plate, we find 

whence , 

u)j «: 6.9(D/moa*) 

This is actually a poor estimate, the correct coefficient being 
approximately 10.22 rather than 6.9. The reason lies essentially in the 
fact that the sum includes a double infinity of terms corresponding to 
both radial and circumferential nodal lines and the frequencies are not 
well separated. A better approximation may be obtained with great labor 
by considering the EME2, or, much more easily, by excluding all modes 
having radial nodal lines. This is accomplished by considering a unit 
load uniformly distributed on the circumference of a circle of radius p. 

The deflection is (1*1) 

z(r,p) = [(r 2 +p 2 )log g (p/a) + (l+p 2 /a 2 ) (a 2 -r 2 )/2 - (p 2 -r 2 ) J/8ttD 
for r £ p, and 

z(r,p) = C(r 2 +p 2 )log fe (r/a) + (l+p 2 /a 2 )(a 2 -r 2 )/2]/8TrD 
for r £ p. The mass function is 
m(r) = 2 , imn 0 

and without difficulty we obtain the DME1 and the DME2 



2 

i=l 



oj7 =(mo/l6TrDa 2 ) 



( 2v f a 

(a >- 



r 2 ) 2 rdrd0 = m 0 a , */*l8D 
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2w 7 2 - moaV96D; 2C * (17/184320) (nioaVD)* 
i=l 1 i=l 



from which we get the respective estimates 

w x » 9 . 8 (D/mo a H )t « t a 10 . 2042 ( D/nioa 4 f L 

It is difficult to assess the accuracy of the latter result. 

Ihe correct coefficient is the square of the smallest positive root 

of the equation 

J 0 (x) Ij(x) + Jj (x) I 0 (x) = 0 

2 

which is given by Rayleigi (9) as (3.20) » 10.24. However, an upper 

bound given by Timoshenko (13) is 

[19.2(51-2(519)* ] 4 « 10.217 



so that our result is definitely in error by no more than 0.13$. 



Example 8 Consider a uniformly tensioned 
uniform circular elastic membrane of radius 
a. The influence function Is 

z(P,Q) = z(r,0;p,4>) = 

- <i/o*t) i°e 6 ( a *t r ' p2 : 2a2roco ° (e ~fb 

^ a 2 [r 2 +p 2 -2rpcos (&-<(>)] 
which Is a rearrangement of a form given by 

Ftektorys (10) with the addition of the membrane 

tension T in the denominator. Ihe DME1 gives 




Notation for analysis 
of circular membrane. 



z(P,P)dS = ~ 

Thus, the DME1 is useless except that it demonstrates the divergence of 
the double sum 



Z Z< n 

m=0 n=l * 
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where J is the nth positive zero of J (x). We have been unable to 
m,n m 

evaluate the integral which gives the DME2. Fbr a successful numerical 
evaluation we do as was done in the preceding exanple , namely confine 
attention to axisyirmetric modes. The deflection at radius x of such a 
circular membrane loaded by unit force uniformly distributed on a 
circle of radius y is 

z(x,y) = (1/2 tjT) log e [a/Max(x,y)] 
and the mass function is 

m(x) = 2Trm 0 x 

There is no difficulty in obtaining the DMEL and the DME2, viz. 

Dg QO 

2)w7 2 = moaV^T = n£a 4 /32T 2 

i=l 1 i=l 1 

As in the evaluation in the preceding example, these sums include 
only frequencies for axisymmetrlc modes. The exact results in the 
present case are 

U i = **0 

Since Jq 2. 1 10 1 l8, the DMA1, which gives the coefficient 2, is in error 
by 17 % and the DMA2, which gives the coefficient (32/* is in error by 
only 1 . 1 %. 

Incidentally, we have here established the interesting and possibly 

useful results „ 

00 00 

Zj" 2 . = 1/1J; 2jn\=l/32 

i=l U<1 i=l u » 1 

Example 9 V/e now consider the lateral vibrations of a simply supported 
rectangular elastic plate having dimensions a and b. V/e assume that the 
ratio p = b/a does not exceed unity. The conpliance is (14) 



00 oo 
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Thus, the DME1 is 



2i w7 2 = (^v,a 4 /D)T , J 1 (iti 2 it 2 + n 2 iT 2 p 2 )' 2 

i*l 1 fr'i i 

m=l n=l 

00 



(m 0 a l> / 1 <Tr l *p 4 D)^ [n 2 Tr 2 p 2 csch 2 (mrp) + nirpcoth(rorp) - 2 ]/n 4 



n=l 

which converges fairly rapidly. The DM22 gives 



CO 



00 00 



2 = (moa^/Tr'D) 2 / / (m 2 + n 2 p 2 )"* 

i=l 1 T 1 

m=l n=l 

CO 

= (m|aV 96 p 8 Tr 8 D 2 ) x , [2n 4 TT 4 p 4 ( csch 4 n7Tp+2csch 2 nTrpcothn7rp) 
n=l 



+ 12n 3 TT 3 p 3 csch 2 nTrpcothnirp 



15nTTfKnTTpcsch 2 n7rp+cothmrp-A8]/n 8 



and this series converges quite rapidly. Thus, for the particular case 
a = b, by taking only one term of the series we get 

« 6.7556*10”^moa 8 /D 2 ;w l «19.6l5(D/m o a 4 ) /;i 
The exactly correct coefficient is 2it 2 = 19.739 so that our estimate is 
quite good. 

For this particular problem, a simple formula (17) gives all of the 
double infinity of natural frequencies, viz. 

m = iT 2 (m 2 /a 2 + n 2 A> 2 ) (D/mo/ 1 

in f n 

and the correctness of the DMS1 and the DM22 can be verified term by term. 
It is merely fortuitous that the exact result for u>, ( =tA, i,i) can be obtained 
in this case by taking only one tern of the double series shown above. 
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